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Abstract 

Mixed-substrate microbial growth is among the most intensely studied systems in molecular micro- 
biology. Several mathematical models have been developed to account for the genetic regulation of such 
systems, especially those resulting in diauxic growth. In this work, we compare the dynamics of three 
such models (Narang, Biotech. Bioeng., 59, 116, 1998; Thattai & Shraiman, Biophys. J, 85, 744, 2003; 
Brandt et al, Water Research, 38, 1004, 2004). We show that these models are dynamically similar — 
the initial motion of the inducible enzymes in all the models is described by Lotka-Volterra equations 
for competing species. In particular, the prediction of diauxic growth corresponds to "extinction" of one 
of the enzymes during the first few hours of growth. The dynamic similarity occurs because in all the 
models, the inducible enzymes possess properties characteristic of competing species: Their synthesis is 
autocatalytic, and they inhibit each other. Despite this dynamic similarity, the models vary with respect 
to the range of dynamics captured. The Brandt et al model captures only the diauxic growth pattern, 
whereas the remaining two models capture both diauxic and non-diauxic growth patterns. The models 
also differ with respect to the mechanisms that generate the mutual inhibition between the enzymes. 
In the Narang model, the mutual inhibition occurs because the enzymes for each substrate enhance the 
dilution of the enzymes for the other substrate. In the Thattai & Shraiman model, the mutual inhibi- 
tion is entirely due to competition for the phosphoryl groups. Elements of all the models appear to be 
necessary for quantitative agreement with data. 

1 Introduction 

When microbial cells are grown in a batch culture containing a mixture of two carbon sources, they often 
exhibit diauxic growth, which is characterized by the appearance two exponential growth phases separated by 
a lag phase called diauxic lag 0|. The most well-known example of this phenomenon is the batch growth of 
E. coli on a mixture of glucose and lactose (Figure^). Early studies by Monod showed that in this case, the 
two exponential growth phases reflect the sequential consumption of glucose and lactose ^^|. Moreover, only 
glucose is consumed in the first exponential growth phase because the synthesis of the peripheral enzymes for 
lactose (the enzymes that catalyze the transport and peripheral catabolism of lactose) is somehow abolished 
in the presence of glucose. During this period of preferential growth on glucose, the peripheral enzymes 
for lactose are diluted to very small levels: 6-7 generations of growth on glucose reduce the enzyme levels 
to ~ 1% of their initial values. Thus, the diauxic lag reflects the time required for the cells to build up 
the peripheral enzymes for lactose to sufficiently high levels. After the diauxic lag, one observes the second 
exponential phase corresponding to consumption of lactose. 

The key to the resolution of the glucose-lactose diauxie is clearly the molecular mechanism by which 
the synthesis of lactose-specific enzymes is abolished in the presence of glucose. The first inroads into this 
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Figure 1: Diauxic growth of E. coli on a mixture of glucose and lactose (from [8]): (a) The optical den- 
sity (□) shows two exponential growth phases separated by a diauxic lag (60 < t < 160 min). The levels 
of /3-galactosidase (O), a peripheral enzyme for lactose, remain low until the beginning of the diauxic lag. 
(b) Evolution of the intracellular cAMP levels during the experiment shown in (a) . The intracellular cAMP 
levels (•) during the first phase of exponential growth on glucose (t < 60 min) are similar to the intracellular 
cAMP levels during the second phase of exponential growth on lactose (t > 160 min). 

problem were made by Monod and coworkers who discovered the mechanism for synthesis (induction) of the 
lactose-specific enzymes in the presence of lactose It was shown that the genes corresponding to the 

peripheral enzymes for lactose are contiguous on the DNA, an arrangement referred to as the lac operon. 
In the absence of lactose, transcription of the lac operon is prevented by a repressor molecule, called the lac 
repressor, which is bound to a specific site on the lac operon. In the presence of lactose, transcription of the 
lac operon is triggered because allolactose, a product of /3-galactosidase, sequesters the repressor from the 
operon, thus liberating it for transcription. 1 

Given this mechanism for induction of the lactose-specific enzymes, it seems plausible to hypothesize 
that the glucose-lactose diauxie occurs because transcription of the lac operon is somehow abolished in the 
presence of glucose. These mechanisms are not fully understood [3]. Until recently, there were two models 
for inhibition of lac transcription in the presence of glucose 

1. cAMP activation: This model postulates that a complex consisting of cyclic AMP (cAMP) and catabo- 
lite repression protein (CRP) must bind to a specific site on the lac operon before it can be transcribed. 
When glucose is added to a culture growing on lactose, the cAMP levels somehow decrease, which re- 
duces the binding of the cAMP-CRP complex to the lac operon, thus inhibiting its transcription rate. 

2. Inducer exclusion: According to this model, enzyme IIA glc , a peripheral enzyme for glucose, is dephos- 
phorylated in the presence of glucose. The dephosphorylated enzyme IIA glc inhibits lactose uptake 
by binding to the lactose permease, the transport enzyme for lactose. This reduces the intracellular 
concentration of allolactose, and hence, the transcription rate of the lac operon. 

Experiments by Aiba and coworkers have shown that the cAMP activation model is not tenable @|. The 
cAMP levels are the same during growth on glucose and lactose (Figure Q]d). Moreover, the lac operon is 
not transcribed even if cAMP is added to a culture growing on glucose and lactose. It is now believed that 
inducer exclusion alone is responsible for inhibiting lac transcription. But in E. coli ML30, the activity 

1 A similar mechanism serves to induce the genes for glucose transport I23L Figure 4]. In the absence of glucose, transcription 
of the ptsG gene, which codes for the transport enzyme, IIBC glc , is inhibited because a repressor called Mlc is bound to a 
regulatory site on the gene. Upon entry of glucose, IIBC glc is dephosphorylated. Dephosphorylated IIBC glc sequesters Mlc 
away from the regulatory site on ptsG, thus liberating the gene for transcription. 
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Figure 2: Nondiauxic growth of E. coli (from |2l|): (a) Simultaneous substrate utilization during batch 
growth of E. coli K12 on a mixture of fumarate (FUM) and pyruvate (PYR). The cell density is denoted 
by c (gdw/L) . This growth pattern is observed with several pairs of organic acids (2j|. (b) Growth pattern 
dependent on the history of the inoculum. When the inoculum is grown on glucose (*), the specific growth 
rate on a mixture of glucose and pyruvate is 0.74 1/hr. When the inoculum is grown on pyruvate (■,♦), 
the specific growth rate on the same mixture is 0.56 1/hr. 

of lactose permease is inhibited no more than ~40% at saturating concentrations of glucose |E Table 2]. 
Likewise, Saier and coworkers, who discovered inducer exclusion in S. typhimurium, found that inducer 
exclusion by glucose inhibits the synthesis of the peripheral enzymes for melibiose, glycerol, maltose, and 
lactose by 10-50% [IE Figures 1-2]. This partial inhibition by inducer exclusion cannot explain the almost 
complete inhibition of the genes for the "less preferred" substrates. 

Although the diauxie has dominated the literature on mixed-substrate growth, there is ample evidence of 
non-diauxic growth patterns. This was already evident from Monod's early studies in which he classified his 
mixed-substrate data into two categories 

EE El- Growth on a particular mixture was called diauxie if the 
growth curve showed the diauxie lag, and normal if it showed no such lag. Yet, the phenomenon of normal 
growth was virtually ignored until recently. In the last few years, several studies have shown that both 
substrates can be consumed simultaneously. Figure shows, for instance, that E. coli consumes fumarate 
and pyruvate simultaneously during batch growth. Egli has summarized all known examples of simultaneous 
substrate utilization in a comprehensive review article He notes that, in general, simultaneous substrate 
utilization is observed when both substrates support low-to-medium specific growth rates, and diauxie growth 
occurs when one of the substrates supports a specific growth rate that is substantially higher than the specific 
growth rate on the other substrate. In addition to simultaneous substrate utilization, there is some evidence 
that the substrate utilization pattern can depend on the history of the inoculum, one example of which is 
shown in Figure 0D (see El 

for other examples). 

The phenomenon of mixed-substrate growth is of fundamental importance in molecular biology as a 
paradigm of the mechanism by which the expression of DNA is controlled. It also has profound implications 
for several large-scale biotechnological processes. 2 This has spurred the development of several mechanistic 
models of mixed-substrate growth. Some of these models are inspired by the detailed, but constantly evolving, 
knowledge of the molecular mechanism for the glucose-lactose diauxie |27l Isfl IHlf . The other models appeal 
to the fact that the phenomenon of diauxie growth is ubiquitous — it has been observed in diverse microbial 
species on many pairs of substitutable substrates (i.e., substrates that satisfy the same nutrient requirements) 

2 The large-scale production of chemicals, such as bioethanol and biopolymers, is economically feasible only if they are derived 
from cheap lignocellulosic feedstocks 0- The pretreatment of these feedstocks yields a mixture of hexoses (primarily, glucose) 
and pentoses (primarily, xylose). The cells that ferment these sugars to useful products typically exhibit diauxie growth with 
preferential consumption of hexoses. 
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including pairs of carbon |Rlfil.ll3T|. nitrogen j^, and phosphorus 0| sources, and even among pairs of electron 
acceptors Thus, it is conceivable there exist some general mechanisms driving the dynamics of mixed- 
substrate growth. These general models, which abstract features common to man y, if not all, mixed-substrate 
systems, include the cybernetic model [13, |2^J and several kinetic models |lL 12(1 . 123 ]. 

The original cybernetic model, which was analyzed in 0], cannot capture nondiauxic growth patterns. 
In this work, we compare the general kinetic models developed by Narang et al [ill, liflj. Brandt et al pj, 
and Thattai & Shraiman j^. Hereafter, we shall refer to these as the N-, B- and T-models, respectively. 
We show that 

1. All the models are similar inasmuch as they exhibit the same general class of dynamics. More pre- 
cisely, the equations describing the initial evolution of the peripheral enzymes are special cases of the 
generalized Lotka-Volterra model for competing species. This similarity arises because all the models 
possess the two defining properties of the Lotka-Volterra model for competing species: Autocatalysis 
(the synthesis of the peripheral enzymes for both substrates is autocatalytic), and mutual inhibition 
(the peripheral enzymes for each substrate inhibit the synthesis of the peripheral enzymes for the other 
substrate). The existence of this similarity implies that the dynamics of the peripheral enzymes are 
analogous to the dynamics of the Lotka-Volterra model. In particular, the prediction of diauxic growth 
by these models corresponds to "extinction" of one of the enzymes. 

2. The models differ with respect to the predicted range of dynamics and the mechanism by which they 
inherit the essential properties of the Lotka-Volterra model. 

(a) The B-model captures only diauxic growth patterns, whereas the N- and T-models capture both 
diauxic non-diauxic growth patterns 

(b) In the N-model, mutual inhibition arises because each enzyme stimulates the dilution of the other 
enzyme. On the other hand, in the T-model, the mutual inhibition occurs because the sugar- 
specific enzymes of the phosphotransferase system compete for phosphoryl groups. 

Comparison with experiments suggests that elements of both all the models are required for capturing the 
data. 



2 The models 

Before describing the models, it is useful to mention a few points. 

1. Although all models contain more or less the same variables, the notation varies considerably from one 
study to another. To facilitate comparison between the models, we have used the same notation for the 
variables. We denote the cells, exogenous substrates, inducers, and peripheral enzymes by C, S%, Xi 
and Ei , respectively. The concentrations of these entities are denoted by the corresponding lower-case 
letters, Sj, e%, Xi and c, respectively. 

2. All the models assume the existence of a small constitutive or background enzyme synthesis rate that 
persists even in the absence of the inducer. We neglect this term since it is generally small compared 
to the induced enzyme synthesis rate. 

3. The cell density and exogenous substrate concentrations are based on the volume of the culture (gdw/L 
and g/L, respectively). In contrast, the concentrations of intracellular variables, such as the enzymes 
and inducers, are based on the dry weight of the cells (g/gdw). 

4. The foregoing choice of units implies that if Z is any intracellular entity produced at the rate, r+ g/gdw- 
hr, and degraded at the rate, r~ g/gdw-hr, then the mass balance for z (in g/gdw) is given by the 
equation 

d(zc) 1 , s dz , (\ dc\ 

— ; — = (rj — r, ) c => — = rj — r. — — — ) z. 
dt v z z ' dt z z \cdtj 

Here, the last term reflects the dilution of Z due to growth. 
We are now ready to describe the key features of the models. 
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Figure 3: Kinetic schemes: (a) N-model (b) B-model [Isj . 

2.1 N-model 

The kinetic scheme for the N-model is shown in Figure Et- It is assumed that 0] 

1. Transport of Si is catalyzed by enzyme E+. The specific uptake rate of Si, denoted r s ^, follows the 
kinetics, r Sti = V s ,ieiSi/(K Si i + si). 

2. Part of the internalized substrate, denoted Xi, is converted to biosynthetic constituents such as amino 
acids and proteins, denoted C~ . The remainder is oxidized to energy (CO2). 

(a) The specific rate of conversion of X, to C~ and CO2 is r X) , = k x< iXi. 

(b) The fraction of Xi converted to C~ is a constant (parameter), denoted Yj. Thus, the specific rate 
of biosynthesis from Xi is Yir x ,i- 

3. The internalized substrate induces the synthesis of Ei. 

(a) The specific synthesis rate of Ei is r &: i = Ve^x™* / (K™\ + a;"*), where m = 1 or 2. 3 

(b) The synthesis of the enzymes occurs at the expense of the biosynthetic constituents, C~ . 

Thus, one obtains the equations 

dsj _ _ = v Si (1} 

j, — r s,iC-, f s i — V s i&i , (I) 

at K s ,i + Si 

77" Ts.i ^xA I 77 1 %i>T'x,i — &x,iXi, (2) 

at \c at ) 

dei ( 1 dc 



dt - \cdtj Cu Te ' 1 ~ Ve < 1 + a* ' (3) 
dc~ , „ / 1 dc\ 

~~dT = ^ + Y ^ ~ ^ + r ^ " U * J c ~ (4) 

These equations implicitly define the specific growth rate and the evolution of the cell density. To see this, 
observe that since all the intracellular concentrations are expressed as mass fractions (g/gdw), their sum 

3 Enzyme induction can be hype rbolic (n, = 1) or sigmoidal (fij = 2), depending on the number of inducer molecules that 
bind to a repressor molecule 
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equals 1, i.e., x\ + x 2 + e\ + e 2 + c = 1, Hence, addition of equations yields 
which can be rewritten in the more familiar form 

2 

= r,.r r,. = 

dt 



^r s ,i-{l-Yi)r x ,i (5) 



r g c, r g 



i=l 

where r g denotes the specific growth rate. 

We can simplify the model by observing that Xi ~ 1CP 3 g/gdw |2|] and r St i,r x> i ~ 1 g/gdw-hr. Thus, 
a;, attains quasisteady state on a time scale of 10~ 3 hr. Moreover, the dilution term r g Xi ~ 10 -3 g/gdw-hr 
is negligibly small compared to r s> i,r X j. Hence, within a few seconds, becomes, « r Sii — r Xji , so that 
r g = J2i r s,i - (1 - ~ X« ^ r s,i- Thus, we arrive at the equations 



(Fir s .i + Y 2 r s . 2 )c, r si = V s<i ei— — (6) 

K- S ,i + Si 

-r s ,iC (7) 



dc 
It 
ds,, 
dt 

de x ni 

= r e ,i - {Yir sA + Y 2 r s ,i)ei, r e>i = V e ,i Rni l + ^ (8) 

V Sli eiSi/(K s<i + 

x t » — - (9) 

c~ = 1 — xi — x 2 — ei — e 2 (10) 

where 10 is obtained by solving the quasisteady state relation, r x ,i ss r s ^ L , for Xi. Substituting (jHJ in the 
expression for r e ^ yields 

which shows that enzyme synthesis is autocatalytic: The larger the enzyme level, the higher its synthesis 
rate. This is a consequence of the cyclic structure associated with the kinetics of induction. Figure^ shows 
that the enzyme, Ei, promotes the synthesis of the inducer, Xi, which in turn stimulates the synthesis of 
even more Ei. This cycle of reactions implies that enzyme synthesis is autocatalytic. 



2.2 B-model 



The B-model is similar to the N-model, the only difference being that the intracellular substrate, Xi, not 
only stimulates the induction of Ei , but also inhibits the induction of Ej (shown in Figure 0Jd as 

arrows with a bar at one end). 4 Assuming that Xi rapidly attains quasisteady state, Brandt et al arrive at 
the equations Q 



dc 
di 
dsj 
~dt 
dej 
dt 



r g ,i + r ga , r 
r, 



'sA^i 's.i 



Yi 



PieiSi/(K si + 



PieiSi/ {K Sy i + sx) + p 2 e 2 s 2 / {K sa + s 2 ) 



(11) 
(12) 
(13) 



4 Brandt et al refer to the intracellular substrate (inducer) and the enzyme induction machinery as signal molecule and 
synthesizing unit (SU), respectively (see Figure 1 of p|). 
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Figure 4: T-model [29] (a) The phosphotransferase system (PTS) for carbohydrate uptake, (b) The kinetic 
scheme. 

where r g i,r St i, Yi denote the specific growth rate, the specific substrate uptake rate, and the yield of biomass 
on the i th substrate, r g denotes the total specific growth rate, and < pi < 1 are parameters called substrate 
preference coefficients. 

Evidently, enzyme synthesis is autocatalytic because of the positive feedback from Xj to Ei. The ap- 
pearance of the specific growth rate, r g , in the expression for r e ^ stems from an additional assumption. It is 
argued that the specific enzyme synthesis rate should be proportional to the specific growth rate to ensure 
that the specific enzyme synthesis rate grows in proportion to the specific growth rate. 

We note finally that Brandt et al scaled the enzyme level, e^, with the level, e*, that would be observed 
during exponential growth on Si alone. Thus, the variable shown in equation ltl3|l corresponds to the 
scaled variable, Ki = e^/e* in (J. 



2.3 T-model 

The T-model is aimed at describing the evolution of the peripheral enzymes belonging to the phosphotrans- 
ferase system (PTS), which catalyzes the uptake of various sugars in certain bacteria The uptake of 
PTS sugars is coupled to their phosphorylation. This is mediated by a cascade of 5 phosphorylation reactions 
involving the successive transfer of a phosphoryl group from PEP to the sugar (Figure QJi) . The first two 
steps, involving phosphorylation of enzyme I (EI) and HPr, are common to all the PTS sugars. The last three 
steps are mediated by sugar-specific enzymes, enzyme IIA (EIIA) and the enzyme IIBC complex (EIIBC), 
which ultimately transfer the phosphoryl group to the sugar during its translocation across the membrane. 
The sugar-specific enzymes are inducible, and their synthesis is coupled since they lie on the same operon. 
Figure QJ} shows the kinetic scheme of the T-model. It is assumed that 

1. There is a maximum flux, J, of phosphoryl groups through the common enzymes, EI and HPr, and the 
sugar-specific enzymes of PTS compete for these phosphoryl group. It turns out that at quasisteady 
state, the specific phosphorylation rate of i th substrate, Ji, (which is equal to the specific substrate 
uptake rate, r Sj i) is given by 

2 

t j T i _ e i S i 

r s,i — Ji — J 1 ) T i — „ i 

1+T1+T2 (H Si + Ci 

where Tj can be interpreted as the demand for phosphoryl groups by the i th sugar — it is an increasing 
function of the exogenous sugar concentration, Si, and the sugar-specific enzyme level, e^. Consistent 
with this interpretation, Ji < J due to competing demands for phosphoryl groups imposed by the 
substrates. 
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2. The specific rate of enzyme synthesis, r e ,i, is proportional to the concentration of the intracellular 
substrate (inducer), Xi, which in turn is proportional to the specific substrate uptake rate, r S) j, Hence 

which implies that enzyme synthesis is autocatalytic. 

3. The substrate concentrations and the specific growth rate are constant — they are treated as param- 
eters. 

Thus, the evolution of the sugar-specific enzymes for the i th substrate, when appropriately scaled, is given 
by the equations 

de l Tj ef Si 

—T = e ii Ti = 7T ■ (14) 

dt 1 + Tl + T 2 Pi Si+ Ci 

Note that the specific growth rate does not appear in the equations, since time is scaled by the parame- 
ter, l/r g . 

3 Results 

We wish to compare the dynamics of the three models described above. At first sight, this seems impossible 
since the T-model describes the dynamics of the enzymes only, whereas the N- and B-models describe the 
dynamics of the enzymes, substrates, and cells. It turns out, however, that the dynamics of the substrates 
and cells are irrelevant on the time scale of interest. Indeed, insofar as the dynamics of mixed-substrate 
growth are concerned, the asymptotic dynamics (f — > oo) of the N- and B-models is of little interest. Much 
more revealing are their dynamics during the first exponential growth phase, since it is these finite-time 
dynamics that determine the substrate utilization pattern. Specifically, diauxic growth will occur if the 
peripheral enzymes for one of the substrates vanishes during the first exponential growth phase. In contrast, 
simultaneous substrate utilization will be observed if the enzymes for both substrates persist during the first 
exponential growth phase. We show below that 

1. In the N-, and B-models, the motion of the enzymes during the first exponential growth phase can 
be described by a reduced system of two equations that are formally similar to the equations of the 
T-model. This makes it possible to compare the N- and B-models with the T-model. 

2. The reduced equations of all the models are different realizations of the generalized Lotka-Volterra 
model for two competing species. Thus, in all the models, the enzymes behave like two competing 
species. In particular, they coexist or become extinct, and these dynamics have meaningful biological 
interpretations in the context of mixed-substrate growth. 

3. The B-model can never capture non-diauxic growth patterns. 

Finally, we compare the mechanisms underlying the dynamics of the N- and T-models. 

3.1 All the models are dynamically similar to the Lotka-Volterra model 

We begin by showing that in the N- and B-models, the dynamics of the enzymes during the first exponential 
growth phase can be described by a reduced system of two equations. To see this, observe that during 
this period, both substrates are in excess, i.e., Si 3> K s ,i- Hence, even though the exogenous substrate 
concentrations are changing, the transport enzymes remain saturated (si/(K s j + Sj) ~1). Now, the cells 
sense the environment through the transport enzymes. Since these enzymes see a quasiconstant environment 
during the first exponential growth phase, they approach quasisteady state levels. It follows that in the N- 
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and B-models, the motion of the enzymes from any initial conditions to the quasisteady state levels can be 
obtained from JHJi an d CSt by replacing Si/(K S; i + Sj) with l. B Thus, we arrive at the reduced equations 

-J7 = V *,i l? \ n, ~ (YiV s . iei + Y 2 V Si2 e 2 ) e» (15) 
at K ei i + e i 

dej 

dt y \piei + p 2 e 2 

Since these equations are formally similar to equation l!14ll . we can compare the dynamics of all three models. 

It turns out that the equations of all three models are dynamical analogs of the generalized Lotka- Volterra 
model for two competing species, which is given by the equations 0, Chapter 12] 

^ = MN u N 2 )N t (17) 

where Ni and fi{Ni 7 N 2 ) denote the population density and specific growth rate of the i th species, respectively, 
and fi{Ni,N 2 ) satisfies the properties 

1. dfi/dN 2l df 2 /dNi < 0, i.e., each species inhibits the growth of the other species. 

2. fi(Ni,N 2 ) < for sufficiently large N\,N 2 > 0, i.e., at sufficiently large population densities, the 
specific growth rates are negative. 

The standard Lotka- Volterra model for competing species is a special case of the generalized model with 

fi{Nx,N 2 ) = n{\ - aaNx - a l2 N 2 ) 

where is the unrestricted specific growth rate of the i th species in the absence of any competition, and 
a.ji, a,i 2 are coefficients that quantify the reduction of the unrestricted specific growth rate due to intra- and 
inter-specific competition [l7j. The analogy between the generalized Lotka- Volterra model and equations 
ill4H16|) becomes evident if we rewrite the latter in the form 

-J - if (ei,e 2 )e 2 , jf ( ei ,e 2 ) = V e ,i ' ni - pW.iei + Y 2 V s , 2 e 2 ) (18) 

d^i rB i \ fB/ ( Pi 



-i = ff { ei ,e 2 )e h ff ( ei ,e a ) = r g ^ 1 (19) 

at \Piei+p 2 e 2 J 

-77- = fi (ei, e 2 )ei, f t (e u e 2 ) = —— e it n = — (20) 

dt 1 + n + t 2 pi Si + e l 

One can check that the functions, , ff , fj satisfy the properties 1 and 2 above. Thus, in all the models, 
the dynamics of the enzymes during the first exponential growth phase are analogous to the dynamics of the 
generalized Lotka- Volterra model for two competing species. 

The dynamics of the generalized Lotka- Volterra model for competing species are well understood 0, 
Chapter 12]. Specifically, the model entertains no limit cycles, so that all solutions ultimately converge to 
some steady state. Despite the absence of limit cycles, the model has a rich spectrum of dynamics. Even in 
the case of the standard Lotka- Volterra model, one can get 4 different types of dynamics depending on the 
parameter values 0|- Indeed, if we define the dimensionless variables, Ui = auN{ and r = r\t, we obtain 
the dimensionless equations 

dUl . _ 012 

— — = (1 - ui - bi 2 u 2 )ui, b 12 = 

dr an 

du 2 . r 2 a 2 i 

—— = p(l ~ b 21 ui - u 2 )u 2 , p = — , b 2 i = . 

dr r% a 22 



5 We have reduced the equations by appealing to intuitive argum ents . This reduction can be justified rigorously by appealing 
to the theorem of continuous dependence on initial conditions (see [2?ij for details). 
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Figure 5: Classification of the global dynamics: (a) The bifurcation diagram of the standard Lotka-Volterra 
model. The full and open circles show stable and unstable steady states, respectively, (b) The bifurcation 
diagram of the B-model. When p\ > P2 (resp., P2 > Pi), -E2 (resp., E\) is rendered extinct during the first 
exponential growth phase. The full and dashed lines show the nullclines for e\ and e2, respectively. The full 
and open circles show stable and unstable steady states, respectively. The arrows show the orientation of 
the vector, {de\j dt.de^j d£), in the regions between the nullclines. 

The steady states of the model are completely determined by the parameters, by, which may be viewed as 
a measure of the extent to which the j th species inhibits the i th species. Figure shows the bifurcation 
diagram of the scaled standard Lotka-Volterra model. The bifurcation diagram shows that when neither 
species inhibits the other species strongly (612, 621 < 1), the two species coexist; when the cross-inhibition is 
asymmetric (612 > 1,621 < 1 or 612 < 1,621 > 1), one of the species is rendered extinct; when both species 
inhibit each other strongly (612,621 > 1), the outcome of the competition depends on the initial condition. 

Given the dynamical analogy between the Lotka-Volterra model and equations l|18H20fl . it is reasonable 
to expect that the peripheral enzymes would yield similar dynamics during the first phase of exponential 
growth. Importantly, these dynamics have simple interpretations in terms of the substrate utilization pattern. 
Indeed, extinction of one of the enzymes during the first phase of exponential growth corresponds to diauxic 
growth; coexistence of the enzymes during this period is the correlate of simultaneous substrate uptake; 
and bistability reflects a substrate utilization pattern which varies depending on the manner in which the 
inoculum has been precultured. 

3.2 The B-model cannot capture non-diauxic growth patterns 

The N- and T-models can capture diauxic growth, simultaneous substrate utilization, and bistable growth 
(see 0, H<| for details). However, the B-model always exhibits diauxic growth. Indeed, the dynamics 
are completely determined by the substrate preference coefficients, p\ and P2 (Figure Eb). If pi > P2, E 2 
becomes extinct during the first exponential growth phase, which corresponds to preferential consumption of 
Si. Conversely, if P2 > Pi, E\ becomes extinct during the first exponential growth phase, which corresponds 
to preferential consumption of S2 ■ 

To see this, it suffices to consider the nullclines of equation ijlfll) . i.e., the curves along which dd/dt = 0. 
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These curves, which separate the eie2-plane into regions in which det/dt is nonzero, are given by the equations 



ei = 0or/f( ei ,e 2 ) = 9 1 = 0, i = 1,2. 

Piei + p 2 e 2 

Evidently, dei/dt — along the e2-axis and the straight line, pid + p 2 e 2 = pi; we shall refer to the latter 
curve as \x. Likewise, de-xjdt — along the ei-axis and the straight line, p\e\ +p 2 e 2 = Pi; we shall refer to the 
latter curve as v. The steady states of Ill9fl lie at the intersection points of the nullclines for e\ and e 2 . Since 
\x and v are parallel, there are no "coexistence" steady states. There are "extinction" steady states at(l,0) 
and (0, l). 6 One can check that the stability of (1,0) and (0, 1) is determined by the disposition of n and 
v. If pi > p 2 , then fi lies above v, and (1,0) is stable, while (0, 1) is unstable. Conversely, if pi < p 2 , then 
v lies above fi, and (0,1) is stable, while (1,0) is unstable. Thus, we conclude that the B-model entertains 
only the diauxic growth pattern. 

3.3 The N- and T-models have different mechanisms of mutual inhibition 

We have shown above that all the models are different realizations of the generalized Lotka-Volterra model 
for two competing species. Furthermore, the B-model cannot capture the non-diauxic growth patterns. In 
what follows, we consider the similarities and differences between the N- and T-models. 

We can develop a better appreciation of the similarities and differences by examining the manner in which 
these models acquire the properties of the generalized Lotka-Volterra model. The latter is characterized by 
two essential properties. 

1. The growth of each species is autocatalytic, i.e., dNi/dt = whenever JVj = 0. 

2. The interaction between the two species is mutually inhibitory, i.e, dfi/dN 2 , df 2 /dN\ < 0. 

It is clear that all the models satisfy the first property precisely because enzyme synthesis is autocatalytic 
(r e ,i = whenever ei = 0). The mechanism that ensures that existence of this property is also identical in 
all the models. It stems from the fact that the enzyme promotes the formation of the internalized substrate 
(inducer) which in turn stimulates the synthesis of even more enzyme. 

The difference between the models lies the mechanism(s) leading to the second property, namely, mutual 
inhibition. In the N-model, each enzyme inhibits the other enzyme by stimulating growth, and thus increases 
the rate of dilution of the other enzyme (e.g., df^ /de 2 < precisely because e 2 appears in the dilution term 
for ei). On the other hand, in the T-model, there is no mutual inhibition due to dilution — in fact, the 
specific growth rate is assumed to be a constant parameter. Instead, each enzyme inhibits the rate of 
synthesis of the other enzyme (e.g., df-f /de 2 < precisely because e 2 appears in the synthesis term for ei), 
and this inhibition occurs due to competition for the phosphoryl groups. 

The N-model has two advantages over the T-model. 

1. It is more general than the T-model since it applies to any pair of inducible substrates, as opposed to 
PTS sugars only. Indeed, the N-model appeals to the two processes — enzyme induction and growth 
— that occur in every system involving inducible substrates. 

2. It explains an important empirical correlation observed in mixed-substrate growth. Based on a compre- 
hensive review of the experimental literature, Harder & Dijkhuizen and Egli have observed that 
in general, both substrates are consumed simultaneously when they support low-to-medium single- 
substrate growth rates. On the other hand, diauxic growth is typically observed when one of the 
substrates supports a much higher specific growth rate. In this case, the substrate supporting the 
higher specific growth rate is usually the "preferred substrate." 

This can be understood in terms of the N-model, wherein each enzyme inhibits the other enzyme by 
enhancing the latter's dilution rate. Thus, enzymes for two substrates that support low-to-medium 
growth rates will coexist since they will not inhibit each other significantly. However, if the two sub- 
strates, say Si and S2, support high and low specific growth rates, respectively, then E\ will strongly 
inhibit the synthesis of E 2 , but E 2 will have little inhibitory effect on synthesis of E\ . Consequently, 
E\ will drive E 2 to "extinction," resulting in preferential utilization of S\. 

6 There is no steady state at (0, 0) since the model is undefined (discontinuous) at this point. 
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The disadvantage of the N-model is that, unlike the T-model, it does not account for inhibition of enzyme 
synthesis. It is conceivable that this occurs by competition for phosphoryl groups. Another mechanism, 
well-documented in the experimental literature, is inducer exclusion |28|. The latter is not accounted for 
by the N- and T-models. Indeed, neither one of these models accounts for direct interaction between the 
enzymes for the two substrates. The effect of the enzymes belonging to the other substrate are exerted 
indirectly by influencing the specific growth rate or demand for the phosphoryl groups. Thus, the N-model 
can be viewed as a general model which is true of every pair of substrates with inducible peripheral enzymes. 
However, for quantitative agreement, it must be modified along the lines of the T-model by accounting for 
specific mechanisms, such as inducer exclusion and competition for phosphoryl groups, that inhibit enzyme 
synthesis. 

It is striking that all the models can predict diauxic growth despite the absence of direct inhibitory interac- 
tions such as inducer exclusion. These dynamics occur precisely because enzyme synthesis is autocatalytic — 
it is this property that makes it feasible for enzymes to become "extinct" during the first exponential growth 
phase. Thus, the models imply that diauxic growth would not be observed if autocatalysis were destroyed. 
This is consistent with the experimental data. Constitutive mutants, in which synthesis of lactose-specific 
enzymes persists even in absence of the inducer concentration (r e ,i| e ._ > 0), do not display the diauxie 0, 
Figure 6]. Similarly, the glucose-lactose diauxie is not observed if the medium contains IPTG, an inducer of 
the lac operon that can enter the cell even in the absence of the lactose permease 0, Figure 7]. 

4 Conclusions 

We compared the similarities and differences between three kinetic models of mixed-substrate growth. We 
showed that 

1. In all three models, the dynamics of the peripheral enzymes are formally similar to the generalized 
Lotka-Volterra model for competing species. This similarity occurs because the peripheral enzymes 
mirror the two essential properties of the Lotka-Volterra model: (a) Synthesis of the peripheral enzymes 
for both substrates are autocatalytic (b) The peripheral enzymes for the two substrates inhibit each 
other. 

2. The model by Brandt et al cannot capture non-diauxic growth patterns. For all parameter values, 
the peripheral enzymes for one of the substrates becomes extinct during the first exponential growth 
phase, thus resulting in diauxic growth. 

3. The models in an d HH capture both diauxic and non-diauxic growth patterns. Both models are 
identical with respect to the mechanism that ensures that peripheral enzyme synthesis is autocatalytic 
— the peripheral enzymes promote the synthesis of the inducer, which in turn stimulates the synthesis 
of even more enzyme. However, they differ with respect to the mechanism that produces mutual 
inhibition. In the Narang model, the mutual inhibition occurs because each enzyme stimulates the 
dilution rate of the other enzyme. In the Thattai & Shraiman model, which applies to PTS sugars 
only, the mutual inhibition stems from competition for phosphoryl groups. 
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